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Dynamical mean field theory combined with finite-temperature exact diagonalization is shown 
to be a suitable method to study local Coulomb correlations in realistic multi-band materials. By 
making use of the sparseness of the impurity Hamiltonian, exact eigenstates can be evaluated for 
significantly larger clusters than in schemes based on full diagonalization. Since finite-size effects 
are greatly reduced this approach allows the study of three-band systems down to very low temper- 
atures, for strong local Coulomb interactions and full Hund exchange. It is also shown that exact 
diagonalization yields smooth subband quasi-particle spectra and self-energies at real frequencies. 
As a first application the correlation induced charge transfer between t^g bands in Nao.3Co02 is 
investigated. For both Hund and Ising exchange the small el Fermi surface hole pockets are found 
to be slightly enlarged compared to the non-interacting limit, in agreement with previous Quantum 
Monte Carlo dynamical mean field calculations for Ising exchange, but in conflict with photoemission 
data. 



I. INTRODUCTION 

Dynamical mean field theory (DMFT)i has been used 
successfully during recent years to describe the electronic 
properties of a variety of strongly correlated materials. 2 
The hallmark of these systems are their complex lat- 
tice geometries, giving rise to intricate single-particle 
properties, accompanied by complex many-electron in- 
teractions in partially filled, nearly localized atomic or- 
bitals. These characteristics lead to a wealth of physi- 
cal phenomena such as metal insulator transitions, ex- 
otic magnetic structures and unconventional supercon- 
ducting phases^ The virtue of DMFT is that it treats 
single-electron and many-electron features on the same 
footing. The key conceptual advance of this approach is 
the mapping of the lattice problem onto an effective im- 
purity problem which is solved numerically exactly. The 
local self-energy is then determined via a self-consistency 
procedure. 

The many-body impurity problem can be solved, in 
the case of single-band systems, by using techniques 
such as numerical renormalization group 4 (NRG), exact 
diagonalization^ (ED), quantum Monte CarloS. (QMC), 
or other schemes^ For computational reasons dynam- 
ical correlations in realistic multi-band materials have 
so far been investigated mainly within QMC 2 , iter- 
ated perturbation theory- (IPT) and the fluctuation ex- 
change method^ (FLEX). The versatility of the QMC ap- 
proach is made feasible by allowing only for Ising-likc ex- 
change interactions to avoid serious sign problems at low 
temperatures* Extensions of QMC including full Hund's 
exchange are presently limited to T = Oifl and rather high 
temperatures (T > 1500 K)iii Spin-flip interactions at 
low finite temperatures can be taken into account in two 
recently developed schemes, namely, the continuous-time 
QMC methodi 2 ^ and the combination of the Hirsch-Fyc 
algorithm with a perturbation series expansion^ 

The aim of this work is to demonstrate that multi- 



band ED/DMFT is a highly useful scheme for the inves- 
tigation of Coulomb correlations in realistic materials. 
In the past, applications of ED/DMFT were limited to 
one-band and two-band systems because of the extremely 
rapid increase of the Hilbert space when the size n s of the 
cluster used to simulate the solid is enlarged>i£ii§*i£ For 
instance, ED with two impurity orbitals, each coupled 
to three bath levels (n s = 8), requires diagonalization of 
matrices with dimension up to 4900 which roughly rep- 
resents the time and storage limit of what is computa- 
tionally meaningful. 18 We show here that a considerable 
simplification of this task can be achieved by exploiting 
the extreme sparseness of the Hamiltonian matrices and 
focusing on the limited number of excited states relevant 
at low temperatures. As a result the cluster size n s = 12 
can be treated at about the same computational cost as 
full diagonalization for n s = 8. This improvement al- 
lows for the first time the application of ED/DMFT to 
realistic three-band systems. 

In this work we study the intercalated layer compound 
Na^CoO^ which exhibits an unusual range of electronic 
properties as a function of Na doping: unconventional 
superconductivity, magnetic and charge ordering, and 
metal insulator transition. Here we focus on one of 
the most puzzling and controversial aspects, namely, the 
topology of the Fermi surface at x = 0.3, the Na con- 
centration at which the system becomes superconducting 
when additionally doped with water^ As a result of the 
octahedral crystal field, the Co 3d bands are split into t^g 
and e g subbands separated by a finite energy gap. With 
Na doping the filling of the t2 g bands can be continuously 
tuned between n = 5 and n = 6. The rhombohedral 
symmetry at Co sites yields a further splitting of ti g or- 
bitals into an a g and two doubly degenerate e' g orbitals. 
According to density functional theory within the local 
density approximation (LDA) the Fermi surface consists 
of a large a g hole surface centered at T and six small e' g 
hole pockets along the TK directionsiSi At x = 0.3 the 0.7 



2 



holes in the t% g band consist of about 0.4 a g holes and 0.3 
e' g holes. Thus, both types of states should be important 
for a variety of electronic and magnetic properties. 

In striking contrast to these LDA predictions, angle- 
resolved photoemission spectroscopy (ARPES) experi- 
ments reveal only the large a g Fermi surface, the e' g 
bands being completely filled. 22 To resolve this discrep- 
ancy several possible explanations have been discussed 
in the literature. Since the onsite Coulomb energy U 
among Co 3<i electrons is about twice the width of the 
tig bands a correlation induced inter-orbital charge trans- 
fer among ti g states could in principle modify the shape 
of the Fermi surface. Dynamical correlations evaluated 
within multi-band QMC/DMFT for a realistic single- 
particle Hamiltonian predict slightly enlarged e' g Fermi 
surface pockets^ The same trend was found in re- 
cent QMC/DMFT calculations^ 2 ^ On the other hand, 
accounting approximately for band narrowing via the 
Gutzwiller approach, filled e' g bands were obtained in 
the strong-coupling U — > oo limit but these calcula- 
tions have not yet been extended to realistic values of U 
and finite exchange integrals J. 

The role of the intercalated Na atoms has also been 
studied. Within LDA, disorder in the Na layer was 
shown to produce potential variations that can localize 
the e' g Fermi surface pockets, at least for large dopingiSi 
A recent LDA+QMC/DMFT study of Na disorder also 
showed that the pockets are likely to disappear at large 
x, but to remain stable near x — 0.3^ 

Although the QMC/DMFT resulto 23 i 24 i 25 so far rep- 
resent the most accurate analysis of correlation induced 
modifications of the Fermi surface of Nao.3Co02, they 
are limited to Ising exchange because of the QMC sign 
problems alluded to above. Since the e' g bands reach only 
about 100 meV above the Fermi level - an energy much 
smaller than typical exchange integrals - there arises 
the question to what extent the inclusion of full Hund's 
exchange might provide an explanation of the observed 
ARPES data. 

To address this issue we have applied ED/DMFT to 
Nao.3Co02- The important result of the present work 
is that dynamical on-site correlations within the Co ti g 
manifold for realistic values of Coulomb and exchange 
energies give rise to a transfer of electronic charge from 
the e'g bands to the a g band. For Hund's exchange and 
T = 10 meV the number of a g holes is decreased to 
0.351 compared to the LDA value 0.40, while the num- 
ber of e'g holes is increased to 0.349 from the LDA value 
0.30. Nearly identical values are found for Ising exchange: 
0.345 (a g ) and 0.355 (e' g ). These results are qualitatively 
consistent with previous QMC/DMFT calculations for 
Ising exchange at T = 30 ... 60 meV. Local Coulomb 
correlations can therefore be ruled out as a possible ex- 
planation of the absence of the e' g hole pockets in the 
ARPES data. 

It remains to be investigated experimentally and the- 
oretically to what extent the photoemission data are 
influenced by surface effects which have played a cru- 



cial role also in other transition metal oxides, such as 
Sr 2 Ru04 and Sr^Cai-^VOs. Also, we note that the 
present ARPES results do not seem to be consistent with 
recent Shubnikov-de Haas measurments of the Fermi sur- 
face of Na . 3 CoO 2 ^ 

The present work establishes multi-band ED/DMFT 
as a useful scheme for describing strongly correlated ma- 
terials. In this sense ED/DMFT can now be regarded 
as complementary to multi-band QMC/DMFT. In com- 
parison to standard Hirsch-Fye QMC treatments, impor- 
tant advantages of ED/DMFT are that spin-flip and pair- 
exchange interactions are fully taken into account, large 
Coulomb energies pose no computational problems, and 
temperatures as low as 5 to 10 meV are readily accessible. 

The outline of this paper is as follows. In Section 
II we present the main ingredients of the multi-band 
ED/DMFT approach. We also argue that it is possi- 
ble within this scheme to evaluate continuous spectral 
functions and self-energies at real frequencies for the ex- 
tended material. Section III presents the application to 
Nao.3Co02 where we focus on the correlation induced 
charge transfer among the partially filled ti g bands. Sum- 
mary and Outlook are given in Section IV. 



II. THEORY 
A. Multi-band ED/DMFT 

Let us consider a material whose single-particle prop- 
erties are characterized by a Hamiltonian HQs). In the 
case of partially occupied ti g bands of a transition metal 
oxide such as Na x Co02, -H"(k) is a 3 x 3 matrix whose ele- 
ments account for direct interactions between ti g orbitals 
as well as indirect interactions via neighboring Oxygen 
ions. In addition we consider on-site Coulomb interac- 
tions which are comparable to or larger than the ti g band 
width. The purpose of single-site DMFT is to derive a 
local self-energy which describes the modification 

of the single-particle bands caused by Coulomb interac- 
tions. The local lattice Green's function is then given by 
the expression 

G a p(iu n ) = (iuj n + (J, - H(k) - T,(iuj n )) (1) 

t — J \ /a/3 
k 

where oj n = (2n + l)ir /ksT are Matsubara frequencies, fj, 
is the chemical potential, and a, f3 = d xy , d xz , d yz denotes 
the tig orbital basis. We consider paramagnetic systems 
and omit the spin index of Green's functions and self- 
energies for convenience. 

In the case of a two-dimensional hexagonal lattice, 
local quantities such as G a p(iu> n ) have only two in- 
dependent elements given by Gn = G22 = G33 and 
Gu = G13 = G23. It is therefore convenient to go over 
to the the a g , e' g basis, where a g = (d xy + d xz + d yz )/V3, 
e ' g i = ( d x V + d xz - 2d yz )/VE and e' g2 = (d xy - d xz )j\[2. 
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Within this basis G becomes diagonal, with elements 

G m= i = G ag = Gn + 2Gi2 (2) 
G m =2,3 = G e ' s — Gn — Gn. (3) 

The reverse transformation is Gn = {G ag + 2G e ' )/3 and 
Gij = {G ag — G e /)/3. Except for the Brillouin Zone 
integral in Eq. we perform all subsequent calcula- 
tions in the diagonal a g , e' g basis denoted by m = 1 ... 3. 
Analogous transformations hold for other lattice struc- 
tures, such as square two-dimensional or cubic three- 
dimensional systems. 

For the purpose of the quantum impurity calculation it 
is necessary to first remove the self-energy from the cen- 
tral site. This step yields the impurity Green's function 



atomic part of the Hamiltonian of the extended mate- 
rial. The intra- and inter-orbital impurity Coulomb en- 
ergies are denoted by U and U'. The exchange integral 
is J. Because of rotational invariance it obeys the rela- 
tion U' = U — 2 J. Spin- flip and pair-exchange terms are 
denoted explicitly by J' . In the case of isotropic Hund 
exchange, one has J' — J. In the case of Ising-like ex- 
change these terms are neglected, so that J' = 0. 

Within the diagonal a g , e' g basis, the cluster Green's 
function has the spectral representationiiSS 



Z ^ E„ 



I(m| C +»I 2 



E„ 



1 J2 [CCrfrn) + G^n)] (8) 



G , m (iw„) = [G m (iu n ) 1 +E m (su n )] x . (4) 

Note that, since G m , Go,™ and E m characterize the ex- 
tended solid, they have smooth, continuous spectra at 
real frequencies. 

Within ED/DMFT^ the lattice impurity Green's func- 
tion Go, TO is approximated via an Anderson impurity 
model for a cluster consisting of impurity levels e m =i...3 
and bath levels Ek = 4 . . .n s coupled via hopping matrix 
elements V m k- Thus, 



Go,m(«w„) « G^ m (iw„) 



where 

Go, m ( iuJ n) = (iUn + /U 



L 

fc=4 



iu n - e k 



(5) 



(G) 



Here, n s = 12 is the cluster size used in the calculations 
discussed below. Since Gq^ is diagonal in orbital indices, 
each impurity level couples to its own bath containing 
three levels. Evidently, in contrast to ImGo, m , at real u> 
Im G$ m is discrete, with the number of poles determined 
by the cluster size. 

The impurity Hamiltonian for the cluster is defined as: 

ma ka 
+ ^2 V ^k[cf na C k a +H.C.] 

aa' J^ma^m' a' 

m m<m' aa' 

— J [ C m~\ C ml C m> l C m' 1 + C m| C m| c )n'T C 'n'l] (^) 

mm' 

where Cml are annihilation (creation) operators for elec- 
trons in impurity level m < 3 with spin a and n ma — 
c mo c mo-, with a similar notation for the bath levels 
k = A...n s . H.c. denotes Hermitian conjugate terms. 
The atomic part of the Hamiltonian is identical with the 



where E v and | v) are eigenvalues and eigenvectors of the 
impurity Hamiltonian and Z = ^2 exp(— f3E u ) is the 
partition function. The excited state Green's functions 
are given by 



G V m A iuJ n) 



MCm»? ^ 

Ey — En — iiU n ' 

\(n\CmaW)\ 2 

^ En — E v - iu n 

li. " 



(9) 
(10) 



In analogy to Eq. 10} the cluster self-energy is given by 
the expression: 

E£(iw n ) = G^ioJn)- 1 - GH^r 1 . (11) 

As will be discussed in more detail in Subsection C, it is 
useful to regard the self-energy as the central quantity de- 
termined via DMFT. The key physical assumption within 
the ED approach is then that the cluster self-energy is 
an adequate representation for the self-energy of the ex- 
tended solid, i.e., 

£ m (i<) wSSlK). (12) 

From Eq. J2J it now follows that 

G m (iu n ) w G m {nu n ). (13) 

The approximate equalities here are important since they 
ensure that, as in the case of Go. m , at real u two alter- 
native representations exist: a continuous one for the 
lattice Green's function G m and self-energy £ TO , and a 
discrete one for the cluster Green's function G^ and self- 
energy EJ^. (See Fig. 13 of Ref. 1 for a nice illustration of 
the fact that continuous and discrete real-w spectra can 
nearly coincide at Matsubara frequencies.) 

Transforming from the a g , e' g basis to the ti g basis 
via £»< = S Qs + 2E e ^, £ i?( y = T, ag - E Cg , we insert this 
self-energy into the solid Green's function Eq. QJ. The 
iteration cycle then has the schematic form: 



G c 



s. 



(14) 



4 



This procedure is repeated until a self-consistent solution 
is found for a given chemical potential. This potential is 
then varied until the calculated total occupancy of the 
tig bands agrees with the desired occupancy. 

In practice, we start from the non-interacting case, 
where the cluster levels e m and Ek and hopping matrix 
elements V m k are found for Go im (iuj n ) — G m {iuJ n ) with 
S m = 0. These cluster parameters are then used to begin 
the iteration at finite Coulomb and exchange energies for 
a given temperature. Finally, the spectral distributions 
are given by 

Am(w) = ImG m (w) . (15) 

The transformation of G m {iuj n ) and S TO (iw n ) to real fre- 
quencies will be discussed in the following subsection. 

At low temperatures, the Boltzmann factors in Eq. JHJ) 
ensure that only a small number of excited states \v) are 
needed. Because of the sparseness of the impurity Hamil- 
tonian matrix, these states can be evaluated exactly by 
using the Arnoldi algorithm^ For a given excited state 
the Green's functions G%^(iu) n ) can then be calculated 
very accurately using the Lanczos procedure, by starting 
from the vectors c^ CT |z/} and c ma \v\ respectively: 

GrnA^n) = p , (16) 

a + - tu n + ai+ _^%» +A „ 

Gmai^n) = p~ ■ (17) 

OQ- - iUn + ai _-J n + bL/ ... 

The reformulation of the cluster Green's function in 
Eq. JSJ), ©, l(TU|) is similar to the one used by Capone et 
ali2£ in the single-band case, except that these authors 
calculate the excited states \v) by applying the Lanczos 
method. This is faster than the exact solution if only 
few excited states are important. However, once exten- 
sive re-orthogonalization is required to obtain sufficiently 
accurate excited states, the Lanczos method becomes in- 
creasingly time-consuming so that the exact approach is 
preferable. 

The Arnoldi algorithm is very useful since it enor- 
mously reduces both storage and time requirements to 
carry out ED/DMFT calculations. For n s = 12 the 
largest Hilbert space for n-f = = 6 has dimension 
N = 853776. Nevertheless, at most M — 23 elements are 
finite in any given row of the Hamiltonian matrix. The 
size of the effective NxM matrix is therefore less than for 
full diagonalization of the largest sector (4900) in the case 
n s = 8. Moreover, at temperatures of about T = 10 meV 
fewer than 30 of the 169 possible (ni-,njj configurations 
have lowest eigenvalues with Boltzmann factors larger 
than 10~ 5 . Most of these sectors contribute only few 
excited states. Only in some (rij,nj) sectors 20 to 40 
eigenstates are important. Thus for a total of about 300 
excited states the Green's functions G^ a (iu! n ) are calcu- 
lated via the Lanczos procedure. Of course, this number 



increases at higher temperatures, and if Boltzmann fac- 
tors smaller than 10 -5 are included for higher precision. 
Finally, the Arnoldi algorithm requires as main compu- 
tational step the matrix times vector operation Hu = v 
where H has effective dimension N x M, the vectors u, v 
are of length N and the internal sum extends only over 
the finite elements of H. Thus, in contrast to full diag- 
onalization, this method is ideally suited for paralleliza- 
tion. Using 16 processors one iteration step takes less 
than 30 Min. 

We have tested the above procedure on the two- 
band model consisting of half-filled subbands of differ- 
ent widths and found excellent agreement with previous 
results for n s = 8 based on full diagonalizationii& 



B. Quasi-particle spectra at real frequencies 

In the ED /DMFT equations given above we explicitly 
distinguish, via the superscript cl, the cluster Green's 
functions and self-energy from the corresponding quan- 
tities of the extended solid. This is done in order to 
emphasize that ED/DMFT is well suited for evaluating 
continuous solid quasi-particle spectra at real frequencies, 
as we now explain. 

The key point is that the solid spectra are to be derived 
from G m {uj), Eq. (JTJ, rather than G^(cj), Eq. ©. Al- 
though in the self-consistent limit both functions agree 
at Matsubara frequencies within some uncertainty de- 
termined by the cluster size, they differ fundamentally 
at real frequencies, with G m (io) exhibiting a smooth 
spectrum and G„(w) consisting of 5 functions. Within 
the spirit of finite-temperature ED/DMFT each iteration 
step involves the evaluation of the solid Green's function 
G m via Eq. Q and the subsequent approximation of the 
associated impurity function Go, m in Eq. Q via the dis- 
crete cluster spectrum of Gg' m , Eq. |JBJ. At a given tem- 
perature, the quality of this projection onto the set of 
Matsubara frequencies depends strongly on the cluster 
size, i.e., on the number of poles at real tu. 

The task at the end of the iteration procedure, of trans- 
forming the lattice subband Green's functions G m (iu) n ) 
from Matsubara frequencies to real u>, is therefore sim- 
ilar to the one faced in QMC/DMFT. The difference is 
that the ED /DMFT data exhibit finite-size effects rather 
than statistical uncertainties. 

Since the ED/DMFT calculations are free of statis- 
tical effects, a convenient extrapolation of G m (iuj n ) to 
real frequencies can be achieved by using, for example, 
the routine rating based on an approximation in terms 
of rational functions. The polynomial basis of this ex- 
trapolation ensures that at real u) the correct continu- 
ous lattice spectrum is found for ImG m (w), rather than 
a discrete spectrum for fixed n s , as obtained by fitting 
Go, m via the cluster impurity Green's function Gq m . We 
use the extrapolation routine ratint in the next section 
to evaluate the subband quasi-particle spectra A m (u>) for 
Nao. 3 Co0 2 . 
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C. Central role of self-energy 

We like to point out that, rather than back- 
transforming the lattice Green's function from imaginary 
to real frequencies, it is more appropriate to perform this 
transformation on the self-energy. In fact, the self-energy 
may be regarded as the central quantity determined by 
the DMFT self-consistent quantum- impurity treatment. 
Since ultimately one is interested in the solid quasi- 
particle Green's function given in Eq. QJ, the single- 
particle features determined by the Hamiltonian -ff(k) 
can evidently be evaluated at real lu. Thus, only the 
real-w variation of the self-energy is needed. As a result, 
k-resolved spectra can be calculated via 

A(k,cj) = --ImTr(u + fj l -H(k)-T,(ufj . (18) 

The usefulness of transforming S from Matsubara to 
real frequencies can be seen at small U. The spectral 
details of G are then almost entirely of single-particle 
character, whereas the self-energy becomes very smooth. 
Since S is a highly convoluted function, with contribu- 
tions stemming from the occupied and unoccupied re- 
gions of the density of states, back-transforming it to 
real lu becomes progressively simpler at small U. In con- 
trast, the back-transformation of G gets more difficult. 
Obviously, since H (k) is known the back-transformation 
of the single-particle part of G can be avoided altogether. 

It is well known also that different cluster properties 
exhibit different convergence with increasing cluster size. 
Single-particle features, such as the local density of states 
at the impurity in a given energy window, converge most 
slowly. Integrated properties, such as the total energy, 
converge more rapidly. Since the phase space for two- 
particle interactions increases extremely fast with clus- 
ter size, many-body properties such as the self-energy 
exhibit even more favorable convergence properties. On 
the other hand, the cluster Green's function involves both 
single-electron and many-electron features. Only the lat- 
ter converge rapidly with cluster size. These different 
behaviors also suggest that it is useful to regard S w E c/ 
as the key quantity derived within ED/DMFT. 

To extrapolate Y\ m (iuj n ) to real frequencies it is first 
necessary to remove the Hartree-Fock limit given by 
(iu) n ) — S TO (iw„ — + oo). The remainder can then 
be transformed to real lu by using the same extrapola- 
tion routine ratint as for the derivation of G m (uj)m^ As 
in the latter case, this routine ensures that the correct 
continuous lattice self-energy is found, rather the dis- 
crete version appropriate for the cluster self-energy. Since 
back-transformation of single-particle features is avoided, 
and quasi-particlc features are calculated from Eq. Q 
at real lu, finer spectral details are retained than when 
G m (iuJn) is transformed to the real axis. Nevertheless, 
in the following section we show that both schemes yield 
consistent spectra. 

Before closing this section we note that the back- 
transformation of the self-energy from Matsubara to real 



frequencies should also be useful within QMC/DMFT. 
By viewing QMC as a tool for evaluating S, the same 
separation of many-body features from single-particle as- 
pects can be achieved if quasi-particle spectra are eval- 
uated from Eq. at real to. Of course, the statistical 
uncertainties associated with QMC make the transfor- 
mation of Yi{iu! n ) more subtle than in the present ED 
case. Nevertheless, after subtraction of the Hartree-Fock 
limit from the real part of S, and after proper nor- 
malization (determined by the asymptotic behavior of 
the imaginary part), the use of the maximum entropy 
method on S(zw„), £(t) should be more accurate than 
on G(iLu n ), G(t), since the back-transformation of single- 
particle features is avoided. 



III. APPLICATION TO NA 3 C00 2 

We have applied the multi-band ED/DMFT approach 
discussed in the previous section to investigate the cor- 
relation induced charge transfer between ti g bands of 
Nao.3Co02- The controversial small Fermi surface pock- 
ets arise from e' g bands that extend less than 100 meV 
above Ep. Since onsite Coulomb energies are much 
larger than the ti g band width, it is important to in- 
quire whether correlations can lead to a filling of the e' g 
bands. 

The single-particle properties for this Na concentra- 
tion were derived from a full-potential linear augmented 
plane wave (LAPW) calculation, which was then fitted 
to a 3 x 3 t2 g tight-binding Hamiltonian including three 
neighboring shells and dda, ddir and ddd matrix elements. 
Further details concerning the LDA band structure are 
given in RefiS. 

For the combined t2 g and e g bands of Nao.3Co02 U 
was calculated to be about 3.7 eV~ For the t^g sub- 
bands a smaller value might be appropriate to account 
for screening via e g electrons. Since precise values of 
U and J are not known, ED/DMFT calculations up to 
U = 5 eV were carried out, assuming J = {7/4. 

Fig. 1 shows the Co a g and e' g subband self-energies 
£ m (zw„) for U = 3 eV and J = 0.75 eV at T = 10 meV. 
Hund exchange is fully included, i.e., J' — J. As a result 
of the large cluster size the spacing between excited states 
of the impurity Hamiltonian is less than 0.001 eV. Al- 
though each impurity orbital hybridizes only with three 
bath levels, the Coulomb and exchange interactions at 
the impurity site induce coupling between all bath lev- 
els, so that the excited states of the total cluster are 
significantly more closely spaced than for n s = 4 in the 
single-band case. Thus, even at temperatures as low as 
T = 2.5 meV finite-size effects are extremely small. The 
slight differences of the real part of the self-energy are 
related to the temperature dependence of the chemical 
potential. 

The a g and e' g quasi-particle weights Z m = 1 / [1 — 
dReS m (u))/du>] ~ 1/[1 — IrnIl m (zu;o)/ci>o] derived from 
these self-energies are Z ag w 0.52 and Z e i w 0.62, giving 
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FIG. 1: Subband self-energies E m (io) n ) of Nao.sCoCh as cal- 
culated within ED/DMFT for U = 3.0 eV, J = 0.75 eV at 
several temperatures: red curves: T = 2.5 meV, blue curves: 
T = 5 meV, magenta curves: T = 10 meV. The Hartree-Fock 
limit is subtracted from ReE m (ioJn). 



effective masses to*^ ~ 1.9 and to*, « 1.6, respectively. 

The corresponding subband lattice Green's functions 
are plotted in Fig. 2 for the same parameters as in 
Fig. 1. As in the case of the self-energies, the distri- 
butions are very smooth and exhibit only very small de- 
viations caused by finite-size effects. 

Fig. 3 shows the variation of the a g and e' g occupa- 
tions with onsite Coulomb energy at T = 10 meV. Both 
isotropic Hund exchange (J' = J) and Ising exchange 
(J 7 = 0) are seen to give nearly identical charge transfer 
from e' g to a g bands. In the present ED/DMFT treat- 
ment this transfer is slightly less pronounced than the one 
found within QMC/DMFT for Ising exchange at higher 
but the trends in both DMFT calculations are con- 
sistent. Thus, both DMFT treatments suggest that local 
Coulomb correlations stabilize the e' g Fermi surface pock- 
ets. The overall topology of the Fermi surface therefore 
remains the same as predicted by the LDA, i.e., correla- 
tion effects cannot explain the absence of the e' g pockets 
from the ARPES data. 

Note that up to U = 5 eV the subband occupations 
show no sign of any reversal of charge transfer. The op- 
posite result found using the Gutzwiller approach^ in 
the U = oo limit might therefore be a consequence of the 
approximate treatment of dynamical correlations. In this 



FIG. 2: Subband Green's functions G m {iu>„) of Nao.3Co02 
for U — 3.0 eV, J = 0.75 eV, at several temperatures: red 
curves: T = 2.5 meV, blue curves: T — 5 meV, magenta 
curves: T = 10 meV. 
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FIG. 3: Subband occupations n ag and n e i (per spin band) 
of Nao.3Co02 as a function of U, with J = U/4 at T = 
10 meV. Solid (red) curves: Hund exchange J' — J, dashed 
(blue) curves: Ising exchange J' = J. The corresponding 
QMC/DMFT results^ for T = 30 ... 60 meV are indicated by 
the symbols (+). The total occupation is 2n ag + 4n e , = 5.3. 



scheme the true complex frequency dependent self-energy 
is replaced by real coefficients representing the shift and 
narrowing of subbands. 

Fig. 4 shows the frequency variation of the subband 
self-energies, as obtained directly from T, m (iu) n ) via the 
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ReS m (w); (b) -ImE ra (u). The Hartree-Fock limit is sub- 
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FIG. 5: Quasi-particle spectra of (a) a g and (b) e' g subbands 
of Nao.3Co02, calculated within ED/DMFT for the same pa- 
rameters as in Fig. 1. Solid (red) curves: spectra derived 
from G m (<jj) via T, m (ui); dashed (blue) curves: spectra ob- 
tained by extrapolating G m (iu;„) to real lu; dotted (black) 
curves: single-particle density of states. 



extrapolation routine ratint. The Kramers-Kronig rela- 
tions between the real and imaginary parts of E m (u;) are 
very well satisfied. The striking asymmetry between the 
negative and positive frequency regions is a consequence 
of the near filling of the subbands. 

As pointed out in the previous section, quasi-particle 
spectra at real frequencies can be derived by back- 
transforming the solid Green's function G m (iuj n ), or by 
first transforming Yi m (iuj n ) and then applying Eq. 
at real oj. The comparison shown in Fig. 4 proves that 
both methods are consistent, and that the latter scheme 
retaines finer spectral details originating from the single- 
particle Hamiltonian. For instance, the e' g spectrum ob- 
tained via Eq. Q and ~E m (oj) shows two peaks below 
Ep which evidently are the shifted and broadened den- 
sity of states features near 0.4 and 0.8 eV below the Fermi 
level. Also, the peak close to Ep exhibits some of the fine 
structure of the single-particle density of states. These 
details are lost if the spectrum is instead derived via back- 
transformation of G m (iu n ). 

Interestingly, the a g spectrum has a steep positive 
slope at Ep, while the e' g spectrum has a strong negative 
slope. This holds true for both methods of evaluating the 
real frequency spectra. Also, the a g quasi-particle spec- 
trum exhibits a pronounced minimum at about 0.1 eV be- 



low Ep. This feature exists already in the single-particle 
density of states and is evidently not obliterated by cor- 
relation effects. 

Both subbands exhibit the usual correlation induced 
band narrowing, quasi-particle damping, and incoher- 
ent spectral weight in the region below the bottom of 
the bands, associated with Hubbard peaks. Qualita- 
tively, these spectra agree with the ones derived previ- 
ously using the QMC/DMFT and the maximum entropy 
method! 23 i 24 i 25 

For comparison with these lattice spectra we show in 
Fig. 5 the cluster spectra obtained by evaluating Eq. © 
at oj + iS instead of ioj n . Because of the large cluster size, 
the level spacing between excited states in the metallic 
phase of Nao.3Co02 is less than 10 -3 eV. Thus, with 
only a minor broadening these cluster spectra look very 
smooth indeed. The spectral weight is located mainly 
in the t^g band region, and there is clear evidence of in- 
coherent weight associated with Hubbard bands. Both 
the a g and e' g spectral distributions exhibit three main 
features in the band region, as in the case of the bare 
density of states. Nevertheless, these features are much 
narrower than in the corresponding lattice spectra shown 
in Fig. 4. This is not surprising since, as pointed out in 
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FIG. 6: Cluster spectra of (a) a g and (b) e' g orbitals of 
Nao.3Co02, as calculated from Eq. JSJ. Solid (red) curves: 
— ImGm((xJ + id) with S — 25 meV; dotted (black) curves: 
single-particle density of states. 



the previous section, the cluster Green's function depends 
on single-particle and many-body aspects. Only the lat- 
ter, represented by the self-energy, converge rapidly with 
cluster size. Thus, increasing the artificial broadening 
of the cluster spectra does not yield the lattice spectra. 
Cluster spectra would resemble more closely those of the 
solid only if n s is significantly increased. To analize pho- 
toemission spectra of the solid material, Eq. should 
therefore be used at real w rather than Eq. JHJ . 

Transforming the subband self-energies S Qs (u) and 
(uj) to the tig basis, they can be used to evaluate the 
momentum dependence of the quasi-particle band struc- 
ture, as indicated in Eq. I|18fl . This is shown in Fig. 7 
for the TK and TM symmetry directions of the hexag- 
onal Brillouin Zone. The e' g Fermi surface pockets are 
associated with the fiat bands reaching just above Ep 
along TK. The maximum of this band is barely shifted, 
but the crossing of Ep occurs slightly farther away from 
the maximum, leading to a weak enlargement of the hole 
pocket. Conversely, the diameter of the large a g pocket 
centered at T is slightly reduced. 

These results demonstrate that, although correlation 
effects induce band narrowing and broadening, they 
do not alter the basic dispersion of the ti g bands of 
Nao.3Co02- In particular, it does not seem possible to 
interprete these bands in terms of a single-band model, 



IV. SUMMARY AND OUTLOOK 

We have applied ED/DMFT to investigate the cor- 
relation induced charge transfer between the ti g bands 
in Nao. 3 Co0 2 . Since previous QMC/DMFT studies of 
this problem were limited to Ising exchange the role of 
full Hund's coupling had remained unresolved. The ED 
calculations show that the charge transfer follows the 
same trend for both Hund and Ising exchange: The e' g 
bands donate some of their charge to the a g bands, i.e., 
the e' g Fermi surface hole pockets are slightly enlarged. 
These ED results fully confirm the trend obtained pre- 
viously within QMC/DMFT for Ising exchange. Thus, 
local Coulomb correlations cannot explain the absence of 
these pockets from the ARPES data. 

These results suggest that ED/DMFT is a useful 
method for the study of strongly correlated three-band 
systems. By exploiting the sparseness of the impurity 
Hamiltonian the cluster size n s = 12 can now be inves- 
tigated, without loss of accuracy, with about the same 
computational effort as n s — 8 using full diagonaliza- 
tion. The large cluster size ensures that the level spacing 
among excited states is very small, so that finite-size ef- 
fects are greatly reduced. 

In view of this improvement, ED/DMFT can now 
be regarded as complementary to QMC/DMFT. We 
point out, however, that in comparison to standard 
QMC/DMFT based on the Hirsch-Fye algorithm, the 
ED/DMFT approach is free of sign problems. Thus, com- 
plete Hund exchange and large on-site Coulomb interac- 
tions can be treated. In addition, very low temperatures 
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can be reached. In principle, results for T — > can also be 
generated within ED/DMFT. This limiting region, how- 
ever, requires special care since certain system properties 
might be sensitive to remaining finite-size effects even for 
n s = 12. 

We have also demonstrated that ED/DMFT provides 
continuous quasi-particle spectra and self-energies of the 
extended solid at real frequencies. Moreover, by con- 
sidering the self-energy as the principal quantity de- 
rived within DMFT, a clear separation of single-electron 
and many-body aspects in spectral distributions can be 
achieved. 

In the future it should be valuable to apply ED/DMFT 
to other multi-band materials that have so far been stud- 
ied only within Ising exchange and QMC/DMFT, but 
where full Hund's coupling is crucial. Also, by exploit- 



ing the fact that the present scheme can easily be paral- 
lelized, much larger cluster sizes are feasible. This opens 
the path for the study of various other material proper- 
ties, including non-local effects. 
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